Post-Newtonian parameters in the tensor-vector-scalar theory 
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We investigate post-Newtonian parameters in the tensor-vector-scalar (TeVeS) theory in a general 
setting while previous researches have been restricted to spherically symmetric cases. Based on the 
assumption that both the physical and Einstein metrics have Minkowski metric at the zeroth order, 
we show 7 = 1 as in the previous researches. We find two remarkable things for other parameters. 
The first is the value @ = 1 while it has been reported that (3^1 for the case when the vector 
field is not purely timelike. This discrepancy occurs from the above assumption which is natural 
as a starting point. The second is the result that the Newtonian potential must be static to be 
consistent with the vector equation. As a result, we cannot determine cvi and 0:2- We consider that 
it is related to the instability against linear perturbation and occurrence of caustic singularities for 
various initial perturbations which have been reported recently. 

PACS numbers: 04.25.Nx, 04.50.Kd, 04.80.Cc 



I. INTRODUCTION 

It is frequently argued that gravitational theories are 
an alternative to dark energy and dark matter. The 
tensor-vector-scalar (TeVeS) theory, which has been pro- 
posed by Bekenstein [IJ , is one of such theories and has 
attracted much attention because of its elegant way to 
explain galaxy rotation curves and the Tully-Fisher law. 
Thus, it is important to investigate its observational con- 
sequences for Te VeS to be a viable theory. 

One of the motivations of TeVeS is to explain the grav- 
itational lensing without dark matter. As for strong 
gravitational lensing, possible explanation using only by 
TeVeS has been reported in [2|. While the difficulty has 
been reported in explaining gravitational lensing of col- 
liding galaxy clusters 1E0657-56 in modified gravities [3|. 
Evolution of the universe and cosmological perturbation 
have also been studied in 0, d, Q . These results tell us 
that TeVeS can explain the observed spectrum of the spa- 
tial distribution of galaxies and of the cosmic microwave 
radiation without dark matter if we include the contri- 
butions by neutrinos and the cosmological constant. The 
scalar field as an alternative to dark energy has been 
pointed out in 0] , if the free function in TeVeS is chosen 
appropriately. 

However, it is not still evident that TeVeS can satisfy 
solar experiments. Parametrized post-Newtonian (PPN) 
parameters has been investigated in [l], Q under the as- 
sumption of spherically symmetric space-time. For this 
reason, although their results show that 7 = 1, it does 
not guarantee the same value in a general setting. As for 
(3, it has been reported that (3 = 1 if the vector field is 
purely timelike and /3 ^ 1 if it is not. It is important to 
notice that the Einstein metric (We explain its meaning 
and the difference from the physical metric in Sec. II.) is 
different from the Minkowski metric at the zeroth order 
in Q . In our opinion, it is the ambiguity peculiar to the 



bimetric theory since the requirement we can safely use 
is that the physical metric is Minkowski metric at the 
zeroth order. However, since the assumption that the 
Einstein metric is also Minkowski metric at the zeroth 
order is simple and hopeful from the aesthetic view, it is 
important to investigate its consequences. It is also im- 
portant to investigate ai and c*2 which characterize the 
preferred frame effects. For these reasons we determine 
PPN parameters based on the general treatment in 
We also follow jToj], which investigate PPN parameters 
in Einstein- Aether theory (ll| . for our calculation and 
presentation. We obtain 7 — (3 = 1 as similar to the 
case where the radial component of the vector field is 
neglected. We cannot determine a\ and a-i because of 
the pathological situation that the Newtonian potential 
must be static. 

This paper is organized as follows. In Sec. II, we ex- 
plain TeVeS and give basic equations for later use. In 
Sec. Ill, we summarize the set up for investigating PPN 
parameters. In Sec. IV, we show the results and explain 
the difference from the previous result Q. In Sec. V, 
consequences and future subjects are discussed. We use 
units in which c = 1 and follow the sign conventions of 
Misner, Thorne, and Wheeler [T^, e.g., (—,+,+,+) for 
metrics. We assume that Greek indices describe space- 
time as a, = 0, 1, 2, 3 and Roman indices describe space 
as i, j = 1, 2, 3. 



II. TEVES AND BASIC EQUATIONS 

The action of TeVeS consists of four parts. We use the 
notation in UM. One is the Einstein-Hilbert action: 



Sg ~ l()7r(V 



-g d^x. 



(2.1) 



We express the metric of this action by <7 M „ (the Einstein 
metric). The second is the vector field's action (gy is a 
dimensionless positive coupling constant): 
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-gd x, 



(2.2) 



which includes a constraint that guarantees the vector 
field to be unit timelike and A is the corresponding La- 
grange multiplier and Uuj := u ff (I — u^^. The third is 
the scalar's action (k is a dimensionless positive param- 
eter while £ is a constant with the dimensions of length, 
and J a dimensionless free function): 



S a 



1 



2k 2 e 2 G 



J (k£ 2 h afS 4)^4>^) V^^x, (2.3) 



where h af3 := g a P — u a u^ with u a := g a ^up. We assume 
that the vector and the scalar fields propagate along the 
Einstein metric. 

One of the points which characterizes TeVeS is the as- 
sumption that usual matter such as electromagnetic fields 
propagate along the physical metric g M „ which is defined 
as 



g a p ■= e 20 g Q/ 3 - 2u a up sinh(20). 



(2.4) 



This point is important to fit the gravitational lensing 
by clusters and galaxies with observations without dark 
matter. Accordingly, the matter action for field variables 
collectively denoted / is: 



Sn 



(W. /£,•••) VHd 4 x 



(2.5) 



where the covariant derivative denoted by | is taken with 
respect to g^ v . 

Variation of the total action with respect to g a @ yields 
the TeVeS Einstein equations for g a p; 



G a 



p = 8ttG \T a p 

+0 a /3 ■ 



(1 



|U(««/3) 



T a f5 



(2.6) 



Here, T^ a up^ := T^up + T^pUa. The sources here are 
the usual matter energy-momentum tensor T a p, the vari- 
ational derivative of S m with respect to g a @ , as well as 
the energy-momentum tensors for the scalar and vector 
fields: 



Taf3 := kG 



1 



3g a 







2k 2 e 2 G 



(2.7) 



3V yg^U^ajU^,^ 

Xu a up, 



g aT g^u^^^.^gap 



where 



p{x) 



dJ(x) 
dx 



(2.8) 



(2.9) 



The equations of motion for the vector and scalar fields 
are, respectively, 



J;/3 



= kG [g aP + (1 + e- 4 ^) u a u ] f a p , (2.10) 

= 8irG (1 - e" 4 *) g^u^f^p . (2.11) 

There is the normalization condition on the vector field: 

u a u a = g aP u a uP = -1. (2.12) 

The Lagrange multiplier A can be calculated from the 
vector equation. 



A = gyUaU^^ -p + ^-u f3 (j)jju'(j)^ 1 
-8vrG (1 - e" 4 ^) vfvPftf . 



(2.13) 



By adjusting 5F, we can explain galaxy rotation curves 
while preserving the Newtonian dynamics near the so- 
lar system. Below, we concentrate on the case which is 
responsible for the Newtonian dynamics, i.e., 



III. SET UP 



(2.14) 



Here, we follow [j| [T(| for our presentation. 

We assume the matter energy-momentum tensor as 



f a/3 = (p + pU + p)v a V p + pg a>3 



(3.1) 



where v a is the four-velocity, p the rest-mass energy den- 
sity, n the internal energy density, and p the isotropic 
pressure of the fluid. 

We assume a nearly globally Lorentzian coordinate sys- 
tem and both the physical metric and the Einstein metric 
are g a p, which is the basic difference from Q, and the 
vector field is purely timelike at the zeroth order. 

As it is usual in the PPN formalism, we consider the 
slow motion and weak gravity limit which are character- 
ized by 

| ~ 0(0.5), p~H~^K) 2 ~0(l) . (3.2) 

We also assume that the components of the physical met- 
ric perturbations h a p := g a p — g a p will be 



'00 



O(l), fc«~0(l), /hk~ 0(1.5). 



(3.3) 



Notice that 0(0.5) terms appear for the consistency with 
JL v l ~ 0(0.5). Notice also that these should be dis- 
tinguished from the Einstein metric perturbations h a p. 
Since we assumed that g a p — r\ a p at the zeroth order, 
we can define h a p := g a p — g a p. This means h a p is same 
order as h a p from the consistency with (|2.4p . i.e., 



0(1) 



/K 3 -~0(1), hoi~ 0(1.5). 



(3.4) 
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This makes important difference from [8j, since we have 
the vector perturbations Su a as 



5u° ~ 0(1), 5u* ~ 0(1.5) , 



(3.5) 



to be consistent with (03). (In [8], 5u l ~ 0(1) is as- 
sumed. One of the most efficient ways to see why Q3.5[) 
appears is to see (14.21) . Lower order of Su a cannot be 
permitted to be consistent with the Lorentz transforma- 
tion and (|3.4|) . I.e., since h oi ~ 0(1.5), u l and u$ are also 
0(1.5). There is a possibility that u l and Ui are lower 
order and {u l — u{) ~ 0(1.5) in some frame. However, 
cancellation in every Lorentz frame is impossible.) As for 
the scalar field, we have 



0-0o + O(l) 



Notice that we can set 



0o = O 



(3.6) 



(3.7) 



which is possible by rescaling the coordinates by factors 
using the relations (|3.3|) and (|3.5[) . This means the 
consistency with the assumption that both the physical 
metric and the Einstein metric are rj a p at the zeroth order 
and characterize the difference from the u r ^ case in 

1- 

The physical metric components are expanded using 
particular potential functions which define the PPN pa- 
rameters as 

.9oo = -1 + 2U - 2/3U 2 - 2£$ w - (Ci - 2£L4 

+(27+2 + a 3 + Ci -2£)$i 

+2(3 7 -2 ( 3+l + C2+£)$2 

+2(1 + C)*3 + 2(37 + 3^-20*4 , (3.8) 



9oj = --( 4 7 + 3 + ai - a 2 + Ci - 2£)Vj 



1 



(1 + a 2 - Ci + 2f)W, 



<? ifc = (1 + 2 7 L%,c ■ 
The potentials are all of the form 

p(y)f 



F(x) = G N J d 3 y 



\x-y\ 



(3.9) 
(3.10) 

(3.11) 



where Gn is the present value of Newton's constant, 
which we determine below in terms of G, k and gy. f 
characterizes each potentials. The correspondences F : f 
are 



U : 1, $i : v lVl , $ 2 : U, $ 3 = H, $ 4 : 

, fa - 2/)j f{y- z ) 



$w : y d 3 zp(z)- 



P 



x-y\ \ \y~ z 



A : 



[«tfa - y) l ] 2 _ w _ u/fa - y)j(x - yf 



(3.12) 



\x~y\ 2 ' ' |a;-y| 2 
The relations useful below are 

F,« = -47rG N/ o/ , (3.13) 



for [/, $1,2,3,4, and V^. 

We also define the superpotential x : 

X = -Gn y d 3 yp|a;- y\ 

which satisfies 

X,H = -2U . 
We also use the relation 

Xfii = Vi-Wi, 



(3.14) 



(3.15) 



(3.16) 



which originally follows from the continuity equation for 
the fluid. 

There is a possibility that goo may depend on x,oo> 
and gij may depend on x,ij- However, it is known that 
these terms can be eliminated using gauge freedoms of 
the coordinate. To obtain this property, we shall impose 
the following gauge conditions: 



h 0iti = AUfi + Bm t i , 



(3.17) 
(3.18) 



where A and B are functions of gy and k which will be 
determined below. 

We also rewrite (12.61) to the form 



R a p — 



8itG (f„ v + (1 - e-^) u~<f l{p u y) + ■ 



(3.19) 



We determine PPN parameters along the following 
recipe: 

(i) Solve the constraint ([2~T2|) for u° to 0(1). 

(ii) Relate with using l|2.4p to 0(1.5). 

(iii) Solve eq.jgH ^ for to 0(1). 

(iv) Calculate eq. (j2~T3f for A t o 0( 1). 

(v) Solve 00-component of eq. (|3.19[) for g on to 0(1). 

(vi) Solve zj-components of eq. (|3.19[) for to 0(1). 

(vii) Solve i-components of eq. (|2.1ip for u l to 0(1.5). 

(viii) Solve Oi-components of eq. (|3.19p forgoi to 0(1.5). 

(ix) Solve eo. fl^TUD for to 0(2). 

(x) Relate /i M „ with h^ v using (|2.4[) to 0(2). 

(xi) Solve 00-component of cq. (|3.19p for goo to 0(2). 
Below, we use the notation as, e.g., 



r _ £00 
l oo — n oo 



U 2 ) A- 



+ (2) , (3.20) 

since we consider quantities to 0(2), where the quantities 
of the superscript (i) mean the quantities of order i. 



IV. CALCULATION OF PPN PARAMETERS 

A. Solving O(l) 
(i): From (|2.12j) . we obtain 

u° = l 



t 00 

2 



(4.1) 



4 



to 0(1) where (|3 .4|) and (13 . 5|1 are used. From this, we 
have 



«o = 9aOU 



-1 



'oo 



(4.2) 



(4.3) 



u% = 9aiU" = h 0i + u l . 

From (|2.12p . we also obtain 

u°- a = , 

to 0(2). 

(ii): From (|2.4j) . we obtain 

h$ =2cf>^ +hg\ h oi = hoi, (4.4) 

using (|2T4]) . f3l6|) . ([3~7]) and (|4"T]) . Below, we also use 

(4.5) 



1,(1) _ _ h oo(i) h (i) _ 

"OO — " ' "ij — 11 



(iii): If we notice that (y/—g), a is the quantity of 0(1) 
or higher, we can evaluate the l.h.s. of (]2 . 10[) to 0(1) as 



1 



-9 



-9(3*0 -uV)0, a ] l/3 



= [(ff o/, -tt a tt^,a],/J = ^ 

where (13.21) are used. Then, we obtain 



(i) 



= kGp 



(4.6) 



(4.7) 



If we compare this with U t a — — AitG^p which is one of 
(|3.13|) . we obtain 



kG 



-U . 



4ttG n 

(iv): From eq. ()2.13p to 0(1), we obtain 



(4.8) 



A = gvu (u 0; ^ - u^) . (4.9) 

Thus, it is necessary to calculate w Q;/3 7 . From ()3.5|) . we 
obtain 



to 0(2). From this general expression, we have 



s : ; _ feotyt 



,(4.10) 



(4.11) 



iH,k + o(^fc,o - hko,i - hi . k ) , (4.12) 



to 0(1.5). Using ||OJ, 631), gTTD and (|4~T3)l to (Oi l. 

we obtain 



\ _ £v«,(i) 

A — 2 00, ii i 



(4.14) 



to 0(1). 

(v): To solve (|3.19p . we should evaluate each terms. 
As for l.h.s. of (|3.19[) , we obtain 



Roo — -^hijhooAj — —hoo,u + (hioj — —ha t o),o 



'^(^00,i) 2 + T^00,j(2/ljj,i — huj) 



to 0(2). Thus, we have 
-Roo = 



2 00,ii ' 



to 0(1). From (j2~7) , we obtain 



(^) 2 

T 00 — „, , TQi = U 



2fcG 



1 

kG 



(4.15) 
(4.16) 

(4.17) 
(4.18) 



to 0(2). We notice that r a( g disappears to 0(1). From 
(|3.1[) . we obtain 

(1 - e-^)^f M(oUo) - 22° (1 - e-^KT M/3 u« 



(4.19) 



to 0(2). Thus, this term also disappears to 0(1). Then, 
if we notice that 



#oo — —A , 9oi — Oij — , 
to 0(1), we obtain 



(4.20) 



(4.21) 



where (pTT]) . P~T4]) and (|3~T5|) are used. Using l[3~3|) and 
(l4~7l). we have 



47T 



(4.22) 



Since the physical metric is g ll u 1 the requirement that we 
recover Newton gravity is expressed as 



h$ := 2U 



(4.23) 



P~2"2"|) and P~2"3"j) tell us that 



to 0(1) and 0(1.5), respectively. Here, we used 
From fO|) . dUTT]), and that is 0(1) or higher 

order quantity, we obtain 



(4.13) 



Gn — G 



2 J 4tt 



(4.24) 



This result is important since the current value of the 
gravitational constant could change its sign depending on 



5 



gv and k. This result is also consistent with [14j except 
the normalization. We also have 



h {1) - 2U 



N 



1 - 



EN 



(4.25) 



(vi): If we use the gauge (13.171) . we obtain Rij = 



-l h Skk- Then > we have 



h (1) 

n ij,kk 



SirGp 



2 J 4tt 



similar to the 00-component. This means 



h {1) -h {1) 5 - h {1) -h {1) 5 - 
— "oo °y i 'Hj ~ "oo °y 



Using (|4.23p . we obtain 



7 



1 



(4.26) 



(4.27) 



(4.28) 



This means that 7 has same value as in the spherically 
symmetric case. Thus, we can conclude that TeVeS has 
not been excluded from the observation related to 7. 



B. Solving 0(1.5) 



(vii): From i-component of eq. (|2.11[) . we obtain 



,i;/3 



u 13 ' 1 .^ to 0(1.5) where we used A ~ 0(1) as 
seen from (14. 14[) . If we remember (|4.3I) and (|4. 13[) . this 
means 



i;0 j:i 
U ,0 = U ,3 



Then, we obtain 



^OoliO + ^Uijj — 2uj 



(4.29) 



(4.30) 



where we used (|4.11l) and (|4. 1 2|) . By taking the spatial 
divergence of this equation, we obtain 



h {1) - 



(4.31) 



Thus, we conclude 

Co - . (4.32) 
This is very important since we obtain 

U, Q = , (4.33) 



from (|4.25| . This means that if we consider the per- 
turbation around Minkowski space-time, the Newtonian 
potential must be static. This does not happen in the 
usual scalar-tensor theory [9( . It is also important to no- 
tice that Einstein-Aether theory has a same property if 
we adopt a Maxwell type vector action |10j. (|4. 33[) also 
means that 



to = X,m = 



(4.34) 



from (|345|) and (|4~8|) . 

(viii): By solving Oi components of eq. (|3.19[) . we ob- 
tain 



hoijj - h j,ij = l&nGpVi 



(4.35) 



where we used ([347)) and (|432|) . Using |348)) . (|430)) and 
(|435|) . we have 

hoj,ij = l'>": ,, ■ (4.36) 
By substituting this to (|4.35|) and using (|4.30p , we obtain 
IGirGpvi — hoijj — Buijj . (4.37) 
Thus, Ui is expressed as 

G 



Bui = h 0l + 4— Vi 



(4.38) 



C. Solving 0(2) 

(ix) : From (|2.4p , we have 

C = 2^) + h® - 2{^f - 2^>fc$ . (4.39) 

This means that we need to solve (12.10)) for (j)^ 2 \ which 
will be performed below, to determine (3. 

(x) : If we notice that 



Too = p(l + n + v lVl - 2U) , 
Tij = pvivj +p5ij , 

r.h.s. of (|2.10[) to 0(2) is evaluated as 



P 



kGp{l + n + 2v lV% ~2U + ^- + 4o - 
If we use (I4.27P , y/—g is evaluated as 



(4.40) 
(4.41) 



») . (4.42) 



l + h, 



(i) 
00 > 



(4.43) 



to 0(1). Using g3J, ((4341) and (14431) . l.h.s. of (12401) to 

0(2) is calculated as 



1 



-.9 



V-9 
+(gV - u l u j )<p :tj 

— <P,n "oo Vii 



- uV>j] + (g ij - uV),i£. 



(4.44) 



Using P~4")) and ([477]) to P~l2]) and ((ITSlj) . we obtain 



,(2) 



kGp(l + 11 + 2viVi -2U + 3 



P 



2h, 



.(4.45) 



(xi): To evaluate r.h.s of eq. (|3.19p to 0(2), we need to 
evaluate each terms to 0(2). Since we have (|4.17l) . (|4. 18f> 



6 



and (|4.19[) to 0(2), the term remained to be evaluated and 
is By looking at (|2.8[) . we notice that we need to 

evaluate A to 0(2). Then, from eq. (|2. 13|1 . we notice that 
we need to evaluate it ^ to 0(2) which is one of the most 
tedious part in our calculation. Using (|4.12[) (quantity of 
0(1.5)) and P~5]) . we obtain 



(4.46) 



From this expression, (|4.12p (quantity of 0(1.5)) and 
we obtain 



0;i 

U .. = ,/ 



0:i 



to 0(2). Using (j4~27)l . we have 



"■00,i , "00 n 00,i 



2 4 
From (|4.46|) and (|4.48|1 . we obtain 



(4.48) 



i;0 1 . . ^ T 



"00 "00,41 



where we used 



1 r 

2 



(4o 3 ) 2 



"00 "oo.m 



(4.49) 



(4.50) 



Then, from ((447)) and (|4"49|) . we can evaluate A to 0(2) 
as 



A = gyuau 1 "'^ - 32TrG(f> (1) p 



-gyUou' .j — 327rG</> ( - 1 

"■00,ij 



Sv 



(4o } ) 2 



"«i,0i ^"OO ft 00,M 



32ttG0 (1) p . (4.51) 



Using this result, we obtain the contribution of 8^ for 
00-component of eq. (|3.19p as 



500 a!i n 

TOO Y Tal19 = 



(4.54) 



Thus, all quantities of r.h.s. of eq. (|3.19p are evaluated. 
As for l.h.s., if we use fOT) . (|4T52"j) and (j4~50|) to gUJ), 

we have 

#oo = ^oo )/l QO,i» _ Tihooji +hio,io - -r (/iqo ) 2 .(4.55) 



Therefore, by summarizing the results ()4.19|) . (|4. 52(1 . 
P~53"l) , P~5l)) and (035}, we obtain 



h-oo.u — 2 
1 

+ 2 



\ YJ ~ ~~2 

(4o ) ) 2 ] ii -2C ) C i -32 7 rG^)p(l-f)" 1 



inGp ( 1 - ^) (1 + n + 2*^ -2U + 3-). (4.56) 



If we use (EL"2"Tj) . ([4351) and (ESI), we obtain 



7(2 



(C) 2 " 

y L 2~ / I"i0,i0 2~ w i,0iJ 

-87rG N p(n + 2v i v l -2U + 3-) 



(0 (1) ) 2 +0 (1) ^ 



2J$h$ M -32nG4>Wp(l 



2 



(4.57) 



If we use (|3.18j) . P~5)) and (|4.25p . we can rewrite this 
equation as 



'oo 
G 



2fcG 
7tGn 



$2-2 — — - C/ 2 + — — 
\47tGn / ttGn 



IT* 



1 - 



2U Z - 8$ 2 - 



2A- 



i--)* 2 



-2 1 



?) >"?)- 

= (2$ 3 + 4$i-4$2 + 6$ 4 ) 



(4.58) 



coo 
ffv 



ffoo 



, 3v c 1 u \ 

+ _ 2 _ l M i,oi — 2 n °°' ii ) ■ 

In a similar way, we obtain 



(C) 2 



+ 16ttG0 (1 V 

(4.52) 



Too - ^-T a pg al3 = ^(foo + fa) 

= ^{l + U + 2viVi -2(7 + 3^) , 
2 p 



(4.53) 



If we remember (I4.38p . we notice that the coefficients of 

~ (2) 

Uj^o should vanish for h Q0 to be expressed without hiQ. 
Thus, we have 



B 



9v 



(4.59) 



Then, we obtain 
2fcG 



kG 



7tGn \47tGn / 7tGn 



u- 



G_ 

Gn 

2$ 3 + 4$! - 4$ 2 + 6$ 4 



77 V 2 



(4.60) 
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As a result, we obtain 



£ = Ci = Ca = C4 = "3 = , 
2 



2/? 



2f- kG 



kG 

7tGn \47tGn 
2 



G_ 

Gn 



1-^ 
2 



(4.61) 



Surprisingly, if we substitute (|4.24[) . we obtain 

(3 = 1. (4.62) 

This is one of the major results in this paper. Thus, we 
cannot restrict TeVeS related to the experiment about (3. 
From the coefficient of $2, we obtain 



2(4-2/? + C2) = -4+^+8(-^ 

7TGn \(jN 

By substituting (|4T24|) . we have 
C2 =0. 



9v 



(4.63) 



(4.64) 



Unfortunately, we cannot determine a\ and a-i since 
we cannot eliminate ui from (|4.38[) . Notice that although 
we apparently obtain 



G 

hoi = -4— U 4 



(4.65) 



for gy = 0, where (|4.4|) is used. We cannot use this 
equation to determine ot\ and since we use (|4.29l) . 
which is obtained from (|2. 1 1 [) under the assumption gy 7^ 
0, to derive (|4^5]) . 



the physical and Einstein metrics are X) a $ at zeroth order, 
it was found that 7 — (3 — 1 as in general relativity. 

We also obtained the important result that the New- 
tonian potential must be static to be consistent with the 
vecvtor equation. This result is pathological to some ex- 
tent. This also happens in Einstein- Aether theory for 
peculiar values of coupling constants [l(|. Because of 
this, we cannot determine a± and a,i. In relation to this, 
it is interesting to remember the instability against lin- 
ear perturbation [HI and occurrence of caustic singular- 
ities for various initial perturbations jl6j |. These results 
also support our pathological situation. Thus, as it has 
been pointed out in these papers [H, [IB] , we should con- 
sider Einstein-Aether type generalization for the vector 
field [HI if we consider TeVeS as an alternative to dark 
matter. It is, of course, very interesting to analyze PPN 
parameters and clarify the conditions for all the PPN 
parameters to be satisfied. This is under investigation. 

Using this generalization, cosmology has been investi- 
gated in [13]. It is also important to investigate neutron 
stars or black holes. In Einstein- Aether theory, these 
features are investigated and show interesting differences 
from the case in general relativity [H, 13 E3] ■ As for 
TeVeS, black holes have been investigated only for the 
cases where the scalar field diverges at the horizon 0, [l3| 
or the scalar field is constant [lq |. There is a possibility 
that the scalar field plays an important role in discussing 
stability of black holes as in black holes with Yang-Mills- 
Higgs fields [2l|. Thus, we also want to consider it as a 
future subject. 
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